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By analyzing B -^ XuH-Ui spectra with a model based on soft-gluon resummation 
and an analytic time-like QCD coupling, we obtain 



O ■ \Vub\ = (3.76 ± 0.13 ± 0.22) x 10" 



where the first and the second error refers to experimental and theoretical errors, re- 
spectively. This model successfully describes the accurate experimental data in beauty 
C^ ' fragmentation, which has similar soft-gluon effects. The \Vub\ value is obtained from 

the available measured semileptonic branching fractions in limited regions of the 
phase-space. The distributions in the lepton energy Eg, the hadron invariant mass 
mx J the light -cone momentum P+ = Ex — \px \ , together with the double distribu- 
tions in {mx,q'^) and {Ei, s™'^^), are used to select the phase-space regions. The q^ is 
the dilepton squared momentum and s™^^ is the maximal m^ at fixed q'^ and E^. The 
\Vub\ value obtained is in complete agreement with the value coming from exclusive B 
decays and from an over-all fit to the Standard Model parameters. We show that the 
slight disagreement (up to +2 a) with respect to previous inclusive measurements is 
not related to different choices for the b (and c) masses, but to a different modelling 
of the threshold (Sudakov) region. 



1 Introduction 

By comparing various spectra in the decays 

5 ^ X„ + £ + z/, (1) 

with the predictions of a model including non-perturbative corrections to soft-gluon dy- 
namics through an effective QCD coupling [1], we obtain a value for the Vub Cabibbo- 
Kobayashi-Maskawa (CKM) matrix element [2, 3] 

|\4^| = (3.76 ± 0.13 ± 0.22) x 10"^ (2) 

where the first and the second error refers to experimental and theoretical errors, respec- 
tively. The model basically involves the insertion, inside standard threshold resummation 
formulae, of an effective QCD coupling (^^(fc^); such coupling is based on an analyticity 
requirement and includes resumming absorptive effects in gluon cascades [4] . A significative 
point is that the model, which has no free unknown parameters, describes rather well also 
5-meson fragmentation data at the Z^ peak [5], where — unlike B decays — accurate 
data are available and there is no uncertainty coming from CKM matrix elements. The 
main properties of the model are sketched in sec. [2l 

We analyze the distributions in the lepton energy Ee, the final hadron invariant mass 
mx, the light-cone momentum P+ = Ex — IpxI, together with the double distributions 
in {mx-i(f) and in (£^£, s™*^"^), with q^ being the dilepton squared momentum and s"^^^ the 
maximal m\ at fixed q^ and Ei. 

The decay rates for the quark-level transitions h -^ uiu^ are proportional to \Vub\'^ 
and two different methods are selected for measuring this matrix element. In the first 
method, one considers a specific exclusive decay by identifying experimentally the final- 
state hadrons. Dynamics is substantially non perturbative and current theoretical predic- 
tions use QCD sum rules, quark models, lattice QCD, etc. The second method involves 
the inclusive hadron final states Xu in eq. ([1]). In general, given a kinematical variable p, 
such as for example the energy Ei of the charged lepton, one measures the number of -B's 
decaying semileptonically to Xu with p in some interval (a, b), divided by the total number 
of produced B's (decaying into any possible final state): 

B[pe a, 6 = — — -— . 3 

IS [B ^ (anythmg)J 

Since the beauty mass mb k. 5 GeV is rather large compared to the hadronic scale, 
rrib ^ A, one can consider semileptonic 6-decays as hard processes, to be treated in per- 
turbative QCD, with inclusive hadron final states coming from gluon radiation. Because 
iKbP/lKibP ~ 10^, 6 — s> cCvi decays constitute a huge background to 6 -^ uiui ones 
as far as inclusive quantities are concernecu. To avoid (or at least substantially reduce) 
such background, one has to consider kinematical regions where 6 — > c transitions are 
kinematically forbidden (or at least strongly disfavored): typically end-point regions. On 
the theoretical side, this restriction has a price, because the available phase-space to QCD 



^ The non-vanishing charm mass reduces the b — > c£iy£ rate roughly by a factor two. 



partons gets strongly reduced. One ends up with the so-called threshold region [6], defined 
as having parametricallja 

mx ^ Ex. (4) 

The perturbative expansion of spectra in the threshold region is affected by large logarithms 
~ ag\o^^{2Ex/mx)^ which must be resummed to all orders in as in order to have a 
reliable result [7-9]. Consistent inclusion of subleading logarithms requires a prescription 
for the QCD coupling in the low-energy region ~ /I — in principle completely arbitrary 
— which in our model is the analyticity condition. Furthermore, Fermi motion, a genuine 
non-perturbative effect related to a small vibration of the h quark in the B meson, comes 
into play when mx becomes as small as ~ \/AEx ■ 

The exclusive determination uses a smaller event sample because it deals with a single 
channel, with the consequence that larger statistical errors are expected. Since the relevant 
hadronic matrix elements can be computed in this case with a first-principle technique, 
namely lattice QCD, one expects that, with increasing computing resources, hadronic un- 
certainties can be systematically (and almost arbitrarily) reduced. On the contrary, the 
inclusive method suffers less from statistics, but needs a modelling of non-perturbative 
QCD effects, which cannot be completely derived from first principles. Asymptotically in 
time, we expect the exclusive measure to take over the inclusive one. 

At present, the determinations from inclusive and exclusive decays are given with a 
relative precision of about 5-8% and 16% [10], respectively, using different computations 
for the inclusive decays and an average of the Lattice QCD determinations for the exclusive 
decays giving 

|Kfe| = (3.51 ±0.21^^6) ^10-3. (5) 

There isa~ +1 — 2cr discrepancy between the inclusive and exclusive measurements, 
depending on the calculation used for the inclusive decays, where the value of \Vub\ obtained 
by the exclusive calculation is always larger than the corresponding measurement obtained 
by the exclusive decays, indicating some "tension" between the above methods. 

A third independent measure of \Vuh\ stems from a general fit of the Standard Model 
(SM) parameters. One assumes the validity of the Standard Model — and therefore also 
the unitarity of the CKM matrix — without using the direct inclusive or exclusive deter- 
minations. The result is [11]: 

\Vub\ = (3.44 ± 0.16) X 10"^ (SMfit). (6) 

The global fit of the SM therefore "prefers" the exclusive determination, while the inclusive 
one is in agreement at ~ 3 a level only. 

It has been suggested that the discrepancy between the value of the experimental mea- 
surement and the inclusive theoretical prediction could signal effects of new physics from 
extra Higgs particles [12]. In our opinion, the above discrepancy does not necessarily imply 
a signal of new physics. In other words, we believe that the "tension" can be dynamically 
explained inside the Standard Model. Even though there are several models in litera- 
ture describing non-perturbative effects in inclusive B decays, which give results perfectly 
consistent with each other [13-18], we believe that a possible interpretation of the above 
scenario is that the theoretical uncertainties have been under-estimated. A re-analysis of 



This region is also called Sudakov region, large-x region and radiation-inhibited region. 



the same data with a rather different model may therefore be usefuL To re-extract \Vub\ 
in this spirit, it is convenient to identify the different dynamical effects which come from 
theory and cannot be extracted from the data. As it will be explicitly shown in sec. E], one 
has to compute, roughly speaking, both: 

1. inclusive rates, (strongly) dependent on the choice of the heavy quark masses rrify (and 
rric), as well as on the QCD coupling at a reference scale (typically asimz)). In sec. |3] 
we present two methods which differ in the treatment of the inclusive quantities; we 
also discuss our choices of the b and c masses; 

2. suppression factors, for the restriction of the kinematical variable p in some exper- 
imentally accessible range. These factors are affected by large threshold logarithms 
and by the related Fermi-motion effects mentioned earlier. 

The discrepancy of our analysis with respect to previous ones does not rely on the estimate 
of inclusive quantities (different choices of quark masses, oi as{mz), etc.), i.e. on point 1., 
but on the modelling of the threshold region, i.e. on point 2. 

In Sec. 3.1 we discuss our choices of quark masses for describing point 1. effects. The 
method used for the point 2. suppression factor is discussed in ref. [1]. Here we just mention 
that within this model the standard perturbative QCD resumming description is modified 
in a minimal way in order to include some infrared non perturbative effects. Since in this 
contest there is no scale separation between perturbative and non perturbative QCD, the 
model is only sensitive to the B meson mass. 

Sec. m contains our results for \Vub\ coming from the various distributions considered 
together with a discussion, while sec. [5] presents our conclusions. 

2 Threshold resummation with an effective coupling 

Let us briefly describe in this section the phenomenological model used to extract \Vub\, 
namely threshold resummation with a time-like QCD coupling for the semi-inclusive B 
decays given in eq. ([T]) (for a more detailed discussion we refer the reader to ref. [1]). 

The first step is the construction of a general analytic QCD coupling from the standard 
one, by means of an analiticity requirement. By requiring that the analytic coupling has 
the same discontinuity of the standard coupling and no other singularity, one obtains for 
example at one loop: 



«5(g^) = ^ 



A' 



(7) 



_iogQ7yi2 g2_^2 

The coupling above has no Landau pole, which has been subtracted by a power correction 
and it is also immediate to check that it has the same discontinuity of the standard one for 
Q^ < 0, i.e. in the time-like region, related to gluon branching. The last term on the r.h.s. 
of eq. ([7]) produces a series of power corrections once it is expanded for Q^ :^ /l^. It is then 
clear that using the effective coupling ([7]) in the standard threshold resummation formula, 
power corrections to the QCD form factor originate from the power corrections in the 
effective coupling through integrations over transverse (fc^) and longitudinal {y) degrees of 
freedom. Since in semi-inclusive decays the gluon is always time-like, we have also included 
the absorptive parts of the gluon polarization function (the well-known "— -ivr" terms) into 



the effective coupling: tliat amounts to a resummation of constant terms to all orders. At 
one-loop for example one obtains: 
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Figure 1: QCD couplings in NNLO. Dotted line (blue): standard coupling asik'jj; contin- 
uous line (red): analytic time-like coupling asi^^). 



Factorization and resummation of threshold logarithms in semileptonic decays leads to 
an expression for the the triple-differential distribution, the most general distribution, of 
the following form [19]: 



d^T 



r dxdwdu 



2E, 



where x = ^^, w = — , u 



C[x,w]as{Q)\(y[u;Q] + D[x,u,w]as{Q)\ 



i-v/i-(2mx/Q)^ 



(9) 



, with El, Ex and mx being the charged lepton 

l+\/l-(2mx/Q)^ 

energy, the total hadron energy and the hadron mass, respectively and the hard scale is 
given hy Q = 2Ex- F = F(as) is the inclusive width of decay ([1]). Furthermore C[x, w; as] 
is a short-distance, process dependent hard factor and D[x,u,w;as] is a short-distance, 
process dependent, remainder function, vanishing in the threshold region u ^ 0. The 
universal QCD form factor for heavy-to-light transitions, resumming to any order in as 
the series of logarithmically enhanced terms to some logarithmic accuracy, has the following 
exponential form in the Mellin moments A^-space [20] go: 



iogMQ)= r^[(i 

Jo y 



y) 



N-l 




yde, 



■Al&siK)] + B[as{Q'y)] + D[as{QV)] 



(10) 



^The Mellin transform of a{u, Q) is as usual aN[Q) = /q (1 — w)^ ^^(u; Q) du. 

''The QCD form factor a[u; Q] has been numerically computed for different values of as{mz) in [1] 



where the functions A{as), B[as) and D[as) have a standard fixed order expansions in as- 
The prescription of our model is simply to replace the standard functions A{as)i B{as) 
and D{as) with the functions A{as), B{as) and D{as) obtained from the standard ones 
by means of the change of renormalization scheme for the coupling constant as ^ 0.3- 

Let us remark that even if our model doesn't contain free parameters to be fitted to 
the data, in a certain sense we have "fitted" the model itself to data. In fact we have 
constructed our model among different possibilities (f.i. different possible prescriptions for 
the low energy behavior of the QCD coupling) (see f.i [8]). A further goal has been to 
describe at the best the very accurate data in B fragmentation [5]. Then, by using the 
discussed perturbative analogy in the resummation formulas, we have made predictions for 
the spectra B decays [1]. 

Threshold suppression — in our opinion the main theoretical ingredient for the measure 
of 1 146 1 5 as discussed in the introduction — is represented by the factor 

W{a,b) = -J- — — =/ p , , , dxdwdu < 1. (11) 

T[B ^ Xui ue\ Ja<p{x,w,u)<b r dxdwdu 

Note that: 

f 1 d^T 

W^(Pmm,Pmax) = / T? , , , dxdwdu = 1, (12) 

Jail i dxdwdu 
where by "all" we mean the whole phase-space of x, w and u. For discussions and results 
on threshold resummed spectra oi B ^ X^ i i^i decays at next-to-leading order see [9] . 

3 1 146 1 extraction: method 

The branching ratio in eq.([3]) can be computed theoretically as: 

B[pe{a,b)] = ^ ^^^-Z ^"^ ^'' "l^ ^^: ^^^ =rsT[B^X^i u,] W{a, h). (13) 
i[B ^ (anythmg)J '- -' 

Since we do not aim at checking QCD but only to extract \Vub\i we can use the experimental 
measure of the B lifetime, which is rather accurate!: 

tb = (1.584 ± 0.007) X 10"^^ s, (14) 

and limit ourself to compute the suppression factor W{a, h) defined in the previous section 
(see eq. ( ITTl) ) and the inclusive B — ^ XuC-Vt rate: 

T[B^X^l v,\ = ^^^^^ F{as) . (15) 

QCD corrections are factorized in the function F{as), which is of the form: 

oo 

F{as) = 1 + 5];F„«^, (16) 

n=l 



^ At the level of accuracy we are interested in, we can safely ignore the difference in lifetimes of the 
neutral and the charged B mesons. One can take for example an average life-time. Similarly, we will take 
average values between charged and neutral B's for the other quantities involved in the rest of the paper. 



with Fn ~ 0{1) numerical coefficientqj. A measure of \Vub\ is provided by direct comparison 
of eq. flT3|) with the corresponding experimental branching ratio. However, within this 
method, a large uncertainty comes from the dependence on the fifth power of the 6-quark 
mass in eq. ( TT5l) l1. 

One can ehminate the above (undesired) dependence on ml, and the related uncertainty, 
by expressing the branching ratio as follows: 

B[pEia,b)] = -^^Wia,b), (17) 

t + l<-c/u 

where we have defined the semileptonic branching ratio {i is a fixed lepton species; in 
practice i = e, n): 



B: 



r [B ^ (anything)] ^ ' 

and the ratio of [b -^ c)/{b —>■ u) semileptonic widths: 

Since B^l is rather well measured, one can give up on its theoretical calculation and replace 
it with the experimental determination [21]: 

BsL = 0.1066 ± 0.0020. (20) 

As discussed in the introduction, W{a, b) strongly depends on the modelling of the thresh- 
old region and is calculated according to the model whose main properties have been 
summarized in sec. |2l for the kinematical distributions listed in sec. HI 

This method presents no m\ dependence, since one has to compute only the ratio of 
widths IZc/u and not the absolute widths. The semileptonic b ^ c width is conveniently 
written as: 

r{B ^XJ u,) = ^^^"^^1^^"'' I{p) F{as) G{as, p) , (21) 

where 

2 

mf. 
p^^^O.l. 22 

mi 

The function I{p) accounts for the suppression of phase-space because of mc 7^ [22] : 

I{p) = l-8p+ 12p2 log - + 8p3 - p^ (23) 

P 

Note that there is an (accidental) strong dependence on the charm mass m^c, because of 
the appearance of a large factor in the leading term in p, namely — 8. As far as inclusive 
quantities are concerned, the largest source of theoretical error comes indeed from the 
uncertainty in p. Most of the dependence is actually on the difference m,h — rric, which can 



^We can safely take m„ = 0. 

''For example, a tiny uncertainty of ±2 % on the rrib mass, which is compatible with the actual estimates, 
corresponds to about a ±10% uncertainty in the total semileptonic b ^ u width, translating into a ±5% 
uncertainty on \Vub\- 



be estimated quite reasonably within the Heavy Quark Effective Theory (HQET) — see 
next section. Finally, the factor G{as, p) contains corrections suppressed by powers of as 
as well by powers of p: 

oo 

G{as,p) = l + J2Gn{p)a^, (24) 

n=l 

with Gn{0) = 0. Note that G(0, p) = G{as, 0) = 1 . By inserting the above expressions for 
the semileptonic rates, one obtains for the perturbative expansion of TZc/u'- 

IT/ 12 



7^c/n = n,/u{p,as,\Vub\/\V,t\) = }^I{p)G{as,p). (25) 



Let us stress that this method actually provides a measure of the ratio |Kjb|/|\4fe|, but 
since the error on \Vcb\ is rather small and theoretically well understood, one is basically 
measuring |Vufc||f| 

We prefer to use the method based on the computation of the ratio of (b ^ c)/{b ^ u) 
semileptonic rates TZc/u (see eq. ( 1T9|) ). instead of the method involving the absolute b -^ 
uiui rate (see eq. ( IT5l) ). because in the former case only ratios of inclusive widths need to 
be evaluated. As a consequence, while eq. (TT5l) depends on the absolute value of m^, eq. (TT9|) 
depends mostly on the difference rUb —rric, opening the possibility of partial cancellation of 
hadronic uncertainties. However, the two methods give compatible values of \Vub\ and the 
method based on eq. fllSp is used to estimate the systematic error on the method based on 
eq. (TT^ . The fact that the two methods give similar results corroborates the conclusion 
that differences with respect to previous analyses originate from a smaller suppression in 
the threshold region in our model compared to previous models, rather than from different 
estimates of inclusive quantities, i.e. different choices of quark masses, QCD coupling, etc. 
In other words, we do not obtain a small value of \Vub\ by an ad-hoc choice of quark masses 
rrib and rric or of asirriz). The difference lies in the different modelling of the threshold 
region. 

3.1 Quark masses 

Since quarks are confined inside observable hadrons, their masses cannot be directly mea- 
sured and their values are biased by the selected theoretical framework. It is well known 
that pole masses are affected by the poor behavior of the perturbative series relating them 
to physical quantities; replacing the pole mass with the MS mass in dimensional regulariza- 
tion slightly improves the asymptotic behavior of the series. Several alternative definitions 
of the 6-quark mass have been introduced in the literature [23] , in order to give a better 
convergence of the first (few) orders of the perturbative series and consequently reduce the 
theoretical errors. 

We have performed the calculation in the MS mass scheme. The MS masses for the 
b and the c quark are taken mbirfib) = 4.20 ± 0.07 GeV and mc{rnc) = 1-25 ± 0.09 GeV 
[21], respectively. However, in order to take into account the uncertainties coming from a 
different scheme definition, we have considered the pole-mass scheme as well. The HQET 



* The average of determinations of | Vcb \ coming from a global fit to the B — > Xclvi and b -^ sj moments 
in the kinetic and IS schemes, in good agreement with each other, is \Vcb\ = (41.6 ± 0.6) x 10^'^ [10,21]. 



in lowest order (static theory) gives the following relation for the difference of the on-shell 
beauty and charm mass: 



rrih — rric = ms — m£, + O 
with 



,., . 1 1 

A 

rUc mb 



(26) 



TUB - ruD ^ 3.41 GeV. (27) 

This value is also consistent with the mass estimates in the kiuEXP and 1 S'exp schemes 
[24]. The leading corrections to ( l26l) . C[yl^(l/?nc — l/mb)], involve the chromomagnetic 
operator and the kinetic one. One can cancel the chromomagnetic corrections by taking 
the "spin-averaged" meson masses [25]: 

ruB + 3 mB*o m/) + 3 m^^o 
rrib - rric^ h ■ ■ ■ , (28) 

where the smaller difference of ~ 70 MeV, implies, for a fixed beauty mass, a heavier charm 

quark masqj: 

ruB + SruB^ _ mD + 3m^.o ^ ^^^ ^^^ 

4 4 

For a consistent subleading computation one also needs the kinetic contribution, which 

cannot be extracted from the data and it is computed theoretically [24,26]. The conclusion 

is that we consider safe to eliminate the c pole charm mass by using a range of rrib — rric 

given by eqs. ([26]) and ([280. 

4 Results 

This is the central section of the paper in which we determine \Vub\ from measured semilep- 
tonic branching fractions, in limited regions of the phase-space, and we perform the corre- 
sponding averages. The experimental analyses are categorized according to the kinematical 
distribution looked at, where selection criteria are applied to define the limited phase-space 
on which the branching ratio is computed. The analyses are: 

1. Ef. where the distribution looked at is the lepton energy (Ei). There are results from 
BABAR [28], Belle [29] and CLEO [30]. The lepton energy ranges from 1.9, 2.0 and 
2.1 GeV for Belle, BABAR, and CLEO, respectively, up to 2.6 GeV. As described 
in ref. [1], we will look only at the range where data are not affected by potential 
b^ciue background: 2.3 GeV < Ee < 2.6 GeV. 

2. mx' where the distribution looked at is the invariant mass of the hadron final state 
[mx)- Both BABAR [31] and Belle [32] have performed one analysis for mx < 1.55 GeV 



^We use niBO = 5279.50 ± 0.33 MeV, tob*o = 5325.1 ± 0.5 MeV, moo = 1864.84 ± 0.17 MeV and 
mD*o = 2006.97 ± 0.19 MeV [21]. 

^"in order to check eqs. (I26p and (|28|) . we have related mt, and nic by considering the ratio between 
B — > Xc T V and B -^ X^ e v branching ratios, which have been measured with a reasonable accuracy at 
LEP and SLD [21]. The corresponding theoretical quantity strongly depends on ttij, and rric [27]: the b and 
c masses extracted in this way corroborates the quark mass relations of the effective theory. 



and mx < 1.7 GeV, respectively. The selection strategy is based on the full recon- 
struction of the B meson on the other side of the event (^-reconstruction analysis). 
The same selected events are looked at also by using two other kinematical distribu- 
tions: (mx, q"^) and P+, as we will see below. In all the three cases a lower cut on the 
photon energy at 1 GeV is applied. 

3. P+: where the distribution looked at is P+ = Ex — Ipxl, Ex and px being the energy 
and the magnitude of the 3-momentum of the hadronic system, respectively. The 
P-reconstruction analysis is used. There are both BABAR [31] and Belle [32] results. 
Both analyses require P+ < 0.66 GeV. 

4. {mx-,q^)'- where the distribution looked at is a two dimensional distribution in the 
plane of the hadronic mass and the transferred squared momentum q^ to the lepton 
pair. The analyses from BABAR and Belle are described in [31] and [32], respectively, 
using the P-reconstruction analysis. Moreover, another technique, called simulated 
annealing, is used by Belle to select events in the {mx,q^) plane [33]. Events with 
hadron mass lower than 1.7 GeV and q^ larger than 8 GeV^ are selected in all the 
cases. 

5. {E^, s^^^): where the distribution looked at is a two dimensional distribution in the 
electron energy and s™^^, the maximal m\ at fixed q^ and E^. There is a result 
from BABAR [34]. The requests on the kinematic variables are E^ > 2.0 GeV and 
sf^^^ < 3.5 GeV^ 

We compute \Vub\ for each of the analyses starting from the corresponding partial 
branching fractions. Then, we determine the average \Vui,\ value using the HFAG method- 
ology [10], together with Ref. [35]. 

Table [1] reports the extracted values of \Vub\ for all the analysis methods and their 
corresponding average. The errors are experimental (i.e. statistical and systematic) and 
theoretical, respectively. The average is: 

\Vub\ = (3.76 ± 0.13 ± 0.22) x 10"^ (30) 

consistent with the measured value of \Vub\ from exclusive decays [10] and the indirect mea- 
sure [11]. The correlation among the analyses has been taken into account when performing 
the average [111. 

The table shows also the criteria used for the determination of the partial branching 
ratio (AB). The \Vub\ values and the corresponding average are plotted in Figure [2j 

Several sources of theoretical errors have been considered: 



• in addition to our preferred method based on eq. fITTj) . we also use the method based 
on eq. flTSj) to extract the value of \Vub\- Since the two methods described in the 
previous section basically involve different inclusive quantities, this error allows a 
cross-check of their evaluations, i.e. basically of the choices of the b and c masses 
adopted; 

"'^"'^ Concerning the B-reconstruction analysis, we use only the mx analysis in Table [U according to the 
approach in Ref. [10]. 
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Table 1: The first column in the table shows the analyses used in the average, the second 
column shows the corresponding values of\Vub\, and finally the last column shows the criteria 
for which AB is available. The final row shows the average value of\Vub\- The errors on the 
\Vub\ values are experimental and theoretical, respectively. The experimental error includes 
both the statistical and systematic errors. 



Analysis 


\Vub\ (10-3) 


AB criteria 


BABAR (Ei) 28 
Belle (Ee) [29] 
CLEO (Ee) [30] 
BABAR (mx) [31] 
Belle (mx) [32] 
BABAR ((mx,g^)) [31] 
Belle imx,q^) [33] 
BABAR {Ee, sjf^^) [34] 
BABAR (P+) [31] 


3.46± 0.14 tof3 
3.25± 0.17 toil 
3.49± 0.20 toil 
4.04± 0.19 tilt 
3.93± 0.26 tiil 
4.14± 0.26 toil 
3.95± 0.42 toil 
3.87± 0.26 toil 
3.45± 0.22 toil 


Ee > 2.3 GeV 

Ee > 2.3 GeV 

Ee > 2.3 GeV 

mx < 1.55 GeV 

mx < 1.7 GeV 

mx < 1.7 GeV, g^ > 8 GeV^ 

mx < 1.7 GeV, g^ > 8 GeV^ 

E^ > 2.0 GeV, sj^^"< 3.5 GeV^ 

P+ < 0.66 GeV 


Average 


3.76± 0.13 toM 





BABAR (E,) 

3.46 ±0.14 ±0.23 

BELLE (E,) 

3.25 ±0.17 + 0.22 -0.21 "" 

CLEO (Ee) 

3.49 ± 0.20 ± 0.23 

BABAR nix 

4.04 ±0.19 ±0.24 

BELLE m^ 

3.93 ± 0.26 ± 0.23 

BABAR m^-q^ 

4.14 ±0.26 ±0.23 

BELLE sim. ann. (m^, q') 

3.95 ± 0.42 ± 0.22 

BABAR (Ee 



niTiaxA 

3.87 ± 0.26 + 0.23 - 0.24 



BABAR P^ 

3.45 ± 0.22 + 0.21 ■ 



0.37 



Average +/- exp + theory - theory 
3.76 ±0.13 ±0.22 



2 



-3i 



IV J [X lO'l 



Figure 2: \Vub\ values for the analyses used in the average and their average. 
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• we compute inclusive quantities both in the MS and pole schemes for the quark 
masses. The \Vub\ value using the MS masses is our default, the value using the pole 
scheme masses, where the ranges defined in sec. 13.11 are used, gives the systematic 
uncertainty. Since in general higher-order corrections are different in the two schemes, 
that should provide an estimate of the size of unknown higher-order effects. 

• we vary the order at which the rate is computed from the exact NLO to the ap- 
proximate NNLO [36]. Since the perturbative series for QCD is believed to be an 
asymptotic one and as = asirrif,) ^ 0.22 in b physics is rather large, that should 
provide a reasonable estimate on the truncation error; 

• we vary all the parameters which enter in the computation of \Vub\ within their errors, 
as given by the PDG [21]. 

What we cannot change is the modelling of the threshold region represented by the factor 
W{a,b), which is fixed in our model, as discussed in the introduction. How good it such 
modelling of the threshold region can only be estimated indirectly, by considering different 
decay spectra, where presumably threshold effects enter in different ways. 

Future work towards improving the determination of the 6-quark mass in the MS and 
the pole scheme will help in reducing our error on \Vub\- 

Table [2] reports the fractional contributions to the theoretical errors due to the different 
sources. A large excursion of the error among the analyses is due to as, which varies from 
a minimum of ±0.6% for the {mx, (f) analyses to a maximum of ±3.5% for the endpoint 
analyses. The largest contributions to the error are given by the charm mass in the MS 
scheme (±4.4%) and the variation from the MS to the pole scheme (from — 1.3% to — 5.2% 
starting from lower to higher masses), due to the conservatively larger error used for the 
pole scheme mass. 

Table [3] shows the \Vuh\ averages for different analysis categories. Note that the {mx, (f) 
analyses tend to have the largest values of \Vub\, while the endpoint analyses the smallest. 

Table 2: The first column of the table shows the different contributions to the theoretical 
errors, the second column shows the corresponding variation, and finally the third column 
shows the percentage contribution with respect to the \Vub\ value. 

Theoretical Errors 



Contribution 


Variation 


Error (%) 


as 


0.1176 ±0.0024 




±0.6 -^ 3.5 


\V,b\ 


(41.6 ±0.6) X 10-3 




±1.4 


rrib ( GeV) 


4.20 ±0.07 




±0.6 


m, ( GeV) 


1.25 ±0.09 




±4.4 


B{B^ XJv,) 


0.1066 ±0.0020 




±1.0 


\Vub\ method 






±0.8 


pole mass ( GeV) 


4.7 < TUb < 5.0, 1.47 <mc < 
3.34 <mb-mc< 3.41 


1.83 


-1.3^ -5.2 


approx. NNLO rate 






±2.0 
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Table 3: The table contains the \Vuh\ values for several analyses and the corresponding 
averages. The errors on the \Vub\ values are experimental and theoretical, respectively. The 
experimental error includes both the statistical and systematic errors. 



Vub\ for endpoint analyses (10 ^) 


BABAR (Ee) [28 
Belle {Ee) [29] 
CLEO (Ee) [30; 


3.46± 0.14 toits 
3.25± 0.17 toil 
3.49± 0.20 till 


Ee > 2.3 GeV 
Ee > 2.3 GeV 
Ee > 2.3 GeV 


Average 


3.43± 0.15 toM 






Vub for "mx analyses (10^^) 


BABAR (mx) [31 
Belle (mx) [32] 


4.04± 0.19 toll 
3.93± 0.26 toil 


mx < 1.55 GeV 
mx < 1.7 GeV 


Average 


4.00± 0.16 toiit 




Vub ioT {mx,q'^) analyses (10^^) 


BABAR {mx,q^) [31 
Belle (mx,g^) [32] 
Belle (mx, q^) [33 


4.14± 0.26 toli 
4.21± 0.37 +[J-i 
3.95± 0.42 toil 


mx < 1.7 GeV, g^ > 8 GeV" 
mx < 1.7 GeV, g^ > 8 GeV^ 
mx < 1.7 GeV, q^ > 8 GeV^ 


Average 


4.13± 0.21 to'tl 






Vub for P+ analyses (10^^) 


BABAR (P+) 31 
Belle (P+) [32] 


3.45± 0.22 toiy 
3.73± 0.32 toil 


P+ < 0.66 
P+ < 0.66 


Average 


3.55± 0.19 t^'ofs 





The larger value of \Vub\ coming from the analysis of the double distribution in (mx, q^) 
is expected on qualitative basis. For mx = 0, the hard scale Q is given by 

q^ 
Q = me - ^^ < 0.71 ruB, (31) 



q 

rriB 



because g^ > 8 GeV . The lower cut on g^ therefore significantly reduces the hard scale 
Q from the "natural" value Q = tub- The point is that our model has been constructed 
to describe P-decay spectra having the (maximal) hard scale Q = m^, and not spectra 
having a smaller hard scale. Indeed, the model was checked against beauty fragmentation 
at the Z^ peak [5] , where the dominant infrared effects are controlled by a hard scale equal 
to niB- In other words, to analyze the {mx-,q^) distribution, we are using the model in a 
region where it has not been checked and it is no surprise that it does not work so well in 
this case. 



5 Conclusions 

We have analyzed semileptonic B decay data in the framework of a model for QCD non- 
perturbative effects based on an effective time-like QCD coupling, free from Landau singu- 
larities. The analysis has considered the kinematical distributions in Ee, mx, and P+, as 
well as the two dimensional distributions in [mx,q'^) and (P^,s™'^^), taking into account 



the experimental kinematical cuts. 
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Our inclusive measure of the \Vub\ CKM matrix element is: 

|\4^| = (3.76 ± 0.13 ± 0.22) x 10"^ (32) 

The errors on the \Vub\ value are experimental and theoretical, respectively. The experi- 
mental error includes both the statistical and systematic errors. 

For the first time, an inclusive value for \Vub\ is obtained which is in complete agreement 
with the exclusive determination. Current literature presents a discrepancy among previous 
inclusive determinations of \Vub\ on one side and the exclusive determinations (~ 2 a") and 
the over-all fit of the Standard Model (~ 3 a) on the other side [10]. 

Let us try to identify the differences between our approach and the previous ones. A first 
difference lies in the selected lepton energy range. According to our analysis, lepton spectra 
below ^ 2.3 GeV measured at the i?-factories suffer from an under-subtracted charm 
background. Because of that, we have limited our analysis to pretty large lepton energies 
Ei > 2.3 GeV. If we take a smaller cutoff, we obtain up to a ~ +2.7% larger average value 
of \Vub\, in order to simulate b ^ c events. As far as theory is concerned, our model seems 
to produce a smaller Sudakov suppression compared to other models constructed on top of 
soft-gluon resummation, such as for example the dressed gluon exponentiation [16]. For a 
fixed experimental rate, a smaller Sudakov suppression implies indeed larger hadronic form 
factors and smaller |Ktfe|'s. As discussed above, our confidence in the model is also based 
on phenomenological grounds; we have checked it in beauty fragmentation [5], where soft 
contributions are similar to those in b decajo. 

In conclusion, the inclusive extraction of \Vub\ requires the calculation of inclusive 
quantities, strongly dependent on b and c masses, as well as the evaluation of threshold- 
suppressed quantities, the latter containing large infrared logarithms and Fermi-motion 
(non-perturbative) effects. We argue that the main difference of our model with respect 
to previous ones is a smaller suppression of the threshold region. 
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